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Abstract 

Accepting validity of self-consistent theory of localization by Vollhardt and 
Woffle, we derive the finite-size scaling procedure used for studies of the 
critical behavior in <i-dimensional case and based on the consideration of 
auxiliary quasi-lD systems. The obtained scaling functions for d = 2 and 
d = 3 are in good agreement with numerical results: it signifies the absence 
of essential contradictions with the Vollhardt and Wolfle theory on the 
level of raw data. The results v — 1.3 — 1.6, usually obtained at d = 3 
for the critical exponent v of the correlation length, are explained by the 
fact that dependence L + Lq with Lq > (L is the transversal size of 
the system) is interpreted as L 1 ^ with v > 1. For dimensions d > 4, 
the modified scaling relations are derived; it demonstrates incorrectness of 
the conventional treatment of data for d — 4 and d = 5, but establishes 
the constructive procedure for such a treatment. Consequences for other 
variants of finite-size scaling are discussed. 



1. INTRODUCTION 

The contemporary situation in investigation of the Anderson localization is character- 
ized by the fact that the results of numerical modelling (see a review article pQ) contradict 
all other information on the critical behavior [TJ [2J [3] . Such situation is unacceptable, since 
undermines a belief in analytical theory. 

The critical behavior of conductivity a and the correlation length £ 

a ocr s , £cx\t\- u (1) 

(r is a distance to the transition point) obtained from the self-consistent theory of local- 
ization by Vollhardt and Wolfle [U [5] , has a form 

, , s = l, 2 < d< oo , (2) 

d > 4 ' ' ' w 
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(d is the dimension of space), and in fact summarizes all known results. Indeed, the formula 
(2): 

(a) distinguishes values d c i = 2 and d c % = 4 as the lower and upper critical dimensions, 
which are known from independent arguments (see [21 [6] for details); 

(b) agrees with theory for d = 2 + e [7J 

v = ^ + 0-e° + 0-e 1 + O{e 2 ); (3) 

(c) satisfies the Wegner scaling relation s = (d — 2)u [8] for d < d C 2] 

(d) gives independent of d critical exponents for d > d C 2, as it is typical for the mean-field 
theory; 

(e) agrees with the results c= 1/2 P [TU] and s = 1 [11J for d = oo; 

(f) agrees with the experimental results s « 1, v « 1 for d = 3, obtained by the 
measurement of conductivity and dielectric susceptibility [I2"| H5].F1 

It is clear that the Vollhardt and Wolfle theory gives at least a very successful approxi- 
mation, satisfying all general principles and reproducing all known results. More than that, 
a suspicion arises that the result (2) is exact [14] H. This conjecture is supported by the 
paper [16] . where Eq. 2 is derived without model approximations on the basis of symmetry 
analysis. 

As for numerical results [UJ-[31J, they can be summarized by the empirical formula 
v ~ 0.8/ (d — 2) + 0.5 [25J, which has the evident fundamental defects. Recent developments 
make a situation even worse giving for d = 3 values v = 1.54 ± 0.08 [24], v = 1.45 ± 0.08 
[2J], v = 1.40 ±0.15 [27], v = 1.57 ±0.02 [22], etc. 

In our opinion, it means the existence of serious defects in the conventional numerical 
algorithms. It is not reasonable to call in question the raw data, which are obtained 
independently by many groups; but it is possible to doubt the algorithms themselves, 
which are not based on any serious theory. In particular, there is a possibility of rough 
violation of scaling [32], or existence of the large characteristic length scale [33| 15]. 

In the present paper, the following approach is accepted. We suppose that the Voll- 
hardt and Wolfle theory (Sec. 2) is correct (there are real grounds for such assumption 
[16]) and derive the quantities which are immediately "measured" in the numerical experi- 
ments. Then comparison can be made on the level of the raw data, avoiding the suspicious 
treatment procedure. 

We restrict the discussion by the popular variant of finite-size scaling based on con- 
sideration of auxiliary quasi-lD systems [34J. Thus, instead of the infinite 3D system one 
consider the system of size L x L x L z , where L z — > oo. Such system is topologically one- 
dimensional and does not possess the long-range order: so the corresponding correlation 

1 The paper [13] is especially interesting, since the experiment is made for the nondegenerate electron 
gas and the influence of interaction can be controlled explicitly. 

2 According to Wegner [IS], the term 0(e 2 ) in (3) is finite and large negative. However, this result 
was derived for the zero-component cr-model, whose correspondence with the initial disordered system is 
approximate and valid for small e; so a difference can arise in a certain order in e. 
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Figure 1: Dependence of the scaling parameter £ld/L on the transversal size L of the 
system. 

length is finite. If £id can be calculated, then its dependence at L — > oo allows to 
registrate phase transitions in the initial 3-D system: it appears, that £id/L — > oo in the 
phase with long-range order and £,id/L — > in the phase with short-range correlations 
[3"H 152] . In the numerical studies, the following scaling relation is usually postulated 



It is based on the assumption that the correlation length £ of the considered d- dimensional 
system is the only essential length scale, so L enters only in the combination If this 
relation is valid, then the quantity £id/D depends on L in accordance with Fig. 1: it remains 
constant at the critical point, while all curves for r > (and correspondingly r < 0) can 
be reduced to one universal curve by the scale transformation. If two curves for t = T\ 
and r = r 2 are calculated, then the scale transformation allows to determine the ratio of 
two correlation lengths. Taking succession T\, r 2 , . . ., one can determine £(Tj) apart from 
numerical factor and to investigate its critical behavior. 

We demonstrate below, that the scaling relation (4) is indeed valid in the limit of large 
£ and L for space dimensions d < 4, while calculation of the scaling function F for d = 2 
and d = 3 shows a good agreement with numerical results (Sec. 3). It signifies that the 
Vollhardt and Wolfle theory is confirmed on the level of raw data. Section 4 clarifies why 




(4) 
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values of the exponent v in the numerical experiments for d = 3 are always greater than 
unity: in the vicinity of the critical point the scaling parameter £id/L behaves as t(L + L ) 
with Lq > 0, which is conventionally interpreted as tL 1 ^ with v > 1. 

For higher dimensions, the scaling relation (4) cannot be correct, and it can be stated on 
the level of a theorem. The problem of the Anderson transition can be exactly reduced to 
the 4 field theory [5J ESI ESI EI] , which is non-renormalizable for d > 4 [38, 39J. Therefore, 
the ultraviolet cut-off (i. e. atomic scale) cannot be excluded from results, and £ is certainly 
not the only relevant length scale. However, it is possible to derive the modified scaling 
relations 

y = F(x) (5) 

with 
and 

Zid ri / r /M-i/3 £ [ln(L/a)] 1/6 

2/ = -T M L M] , a = - 1/2 , d = 4, (7) 

demonstrating incorrectness of the conventional treatment of data for d = 4 and d = 5 [TJ, 
but establishing the constructive procedure for such a treatment (Sec. 5). The modified 
scaling (5) with 



(L/a) 



L [(e/a)^l] 1/2 (Va)" e/2 ' 



can be derived also for d = 4 — e dimensions. It can be used for alternative treatment 
at d — 3, in order to estimate the systematic errors related with possible existence of the 
large length scale. Finally, in Sec. 6 we discuss some consequences of the present analysis 
for other variants of finite-size scaling. 

2. VOLLHARDT AND WOLFLE THEORY 

The Vollhardt and Wolfle theory is based on existence of the diffusion pole in the 
irreducible four-leg vertex Z7jdc'(q), 

— + i/(a;, k + k'J(k + k'j z 

entering the Bethe-Salpiter equation and playing the role of the scattering probability 
Wkk' in the quantum kinetic equation. Neglecting the spatial dispersion of the diffusion 
coefficientH and using the estimate in the spirit of r-approximation, D oc (U)" 1 , where (...) 



3 Such possibility was justified in [16 . Attempts to relate the spatial dispersion of the diffusion coefficient 
with multifractality of wave functions [41] ignore the complex- valuedness of the diffusion coefficient and its 
complicated rearrangement near transition |42) . 
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is averaging over momenta, we come to the self-consistency equation of the Vollhardt and 
Wolfle theory 

d d q 1_J 



D 



U + F J 



(10) 



-iu + D(lo, q)q 2 

It can be obtained by approximate solution of the Bethe-Salpiter equation [4J, or by the 
accurate analysis of spectral properties of the quantum collision operator [16J. It can be 
written in the physically clear form, if coefficients are estimated for weak disorder (which 
is actual for lower dimensions) and a situation near the band center in the Anderson model 
is implied: 

A 2-d f dd Q 1 
W 2 ~ D min + J (2tt)<* (-iu/D) + q 2 ' 

|?|<A 

Here E is the energy of the bandwidth order, W is the amplitude of disorder, A is the 
ultraviolet cut-off, D min is a characteristic scale of the diffusion coefficient corresponding 
to the Mott minimal conductivity. Generally, some monotonic function of W is appearing 
in the left-hand side, but it is not essential for subsequent considerations. 
Let introduce the basic integral 



(11) 



I(m) 



1 



d d q 

(2ix) d m 2 + q 2 



(12) 



which can be estimated for m<Aas 



I(m) 



,d-2 



where 



Cd 



Cd/m 2 d , 
C2 ln(A/m) 
1(0) - c d m a 
1(0) - c A m 2 ln(A/m) 
1(0) - c d m 2 A d - i , 

nK d /(2sm(rrd/2)), 

1/27T, 

7rK d /\2sm(7cd/2)\ 
V(8vr 2 ), 
K d /(d-A), 



d<2 
d = 2 
2<d<4 
d = A 
d>4 

d<2 
d = 2 
2<d<4 , 
d = A 
d>4 



(13) 



(14) 



and Kd = 2 d l n d l 2 T(d/2) is the surface of the (i-dimensional unit sphere divided by 
(2ii) d . The metallic phase is possible, when value of 1(0) is finite, i. e. for d > 2. Accepting 
D = const > for uj — > and specifying r as a distance to transition, one has 



D = D r 



T . 



T 



F 2 

W~ ms 



2-d 



(15) 
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i.e. the exponent of conductivity is unity, in agreement with (2). In the dielectric phase 
we make substitution 

D = -iui\ i = m-\ (16) 
where £ is the correlation length. Then Eq. 11 gives 

E 2 

^ a W^ d = 1 ' 

£~aexp^27r^, d = 2, (17) 

£~a\r\- u , d>2, 

with the exponent v defined by Eq. 2. In what follows, we accept a = A -1 , so a is the 
atomic length scale, not necessary coinciding with the lattice spacing. 



3. SCALING FUNCTIONS FOR D < 4 

3.1. Definition of scaling functions 

For description of quasi-lD systems it is sufficient to present the basic integral (12) in 
the following form 

1 ^ r dq\\ 1 



where d- dimensional vector q — (q 1 , q 2 , . . . , qa) is replaced by its transversal and longitudinal 
components 

Q± = (5i,52,---,5d-i) , Q\\ = Qd, (19) 

and the first is considered as discrete, running the usual allowed values. The term with 
q± = has divergency m _1 for m — > 0, so the system is always localized. 
After integration over q\\, the following decomposition is convenient: 

T, S 11 A 

1 ( m ) = Td~i arctg— + 

L a 1 Tim m 



1 „ / 1 A 1 A , 

— 7TT 1^ i arctg = - 1 — r arctg: — r + 

- 1 ^ 2 + ?i \/ m2 + 2 



7rL d - 

kxl<A 



+ FTT £ A arct S A = J i( m ) + J 2( m ) + J 3(°) • ( 2 °) 
7T.L" 1 5j_ | q±\ 

\q±\<A 
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We separated the term with q± = 0, while the remaining sum was rearranged by subtraction 
and addition of the analogous sum with m — 0. In the first term one has trivially 

''< m > = lM2^ + (l)}' < 2o > 

The second term can be transformed by taking the limit A — > oo and substituting q± = 
2ns/ L, where s — (s±, . . . , s^-i) is a vector with integer components = 0, ±1, ±2, . . .: 



L d 

™ = *%hwki*r-^- (21) 



The third term can be estimated at L — > oo by the replacement of summation by integra- 
tion. For finite L and d > 2 it has a structure 

/s(0) = A d - 2 j&a + 61 (f )^ + & 2 (I )^ + •••}. (22) 
Substituting (20-22) in the self-consistency equation (11), we have for d > 2 



d-2 



(-) ~ [r + 0(mV)] + O (£) = 6: + ff„(mL) + ^ (23) 
where we replaced 

r=~h* (24) 

in agreement with definition (15), since 60 corresponds to 1(0), calculated in the integral 
approximation. Expressing r through the correlation length £ of the d- dimensional system 
~ \t\ = ±t) and omitting the terms dissapearing at a — > 0, we have 



i.e. the scaling relation (4) between variables £id/L and consisting of two branches. 
For d = 2 we have instead (22) 

/ 3 (0)==^-ln- + 6i + ... (22') 

Z7T a 

and, using the result from Sec. 2 

E 2 1 , £ 




Figure 2: Function H(z) for b\ = in two and three dimensions, 
obtain the scaling relation in the form 

with the previous definition of H(z). The functions H(z) for d = 2 and d = 3 are presented 
in Fig. 2, where &i = was accepted. 

3.2. Two-dimensional case 

For d = 2, the constant &i can be eliminated by the change of the scale for £ (see below) 
and we can accept b\ = 0. The asymptotics of H(z) for z <C 1 is determined by the last 
term in Eq. 26, while for z 3> 1 the sum in Eq. 26 can be replaced by the integral 



H(z) = 



1 

1 2 ^ , (28) 

— — In z + const , z >> 1 

27T 



so we have in variables y = £id/L and x = £/L 



[ (lA)lnx, x> 1 

j const • x , i < 1 
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The relation between x and y for their arbitrary values can be found by the numerical 
calculation of the sum in (26). 

The definition of £i£> and £ in the Vollhardt and Woffle theory does not coincide with one 
used in numerical experiments. In the former case, £ 2 (and analogously £f D ) is defined as an 
average (r 2 ) for the localized eigenfunction ip(r) |16j . In the latter case, one has in mind the 
definition through the asymptotic behavior exp{— r/£} of the correlation functions, since 
£i£> is calculated as inverse to the minimal Lyapunov exponent^; the scale of £ in numerical 
experiments is arbitrary from the very beginning. Therefore, in comparison of theory with 
numerical results the scales of £i£> and £ should be chosen from the best agreement; in the 
log-log coordinates, such fitting reduces to parallel shifts along two axes. The general form 
of the scaling curve is determined without adjustable parameters. 

In Fig. 3, the calculated dependence of Cid/L on £/L is compared with the pioneer re- 
sults by MacKinnon and Kramer [TB] and the subsequent paper by Schreiber and Ottomeier 
j!9j . which is cited as the most detailed investigation of the 2D systems in the framework 
of the given algorithm. 



3.3. Three-dimensional case 

The given definition of the sum 1^(0) implies the choice of cut-off in the form of the 
cilindrical domain (\q±\ < A, \q\\\ < A). It can be also defined for the spherical (\q\ < A) 
and cubical (\qi\ < A) regions: 



|si|<AL/27r 



(sph) 




1 / AL \ 

-arctg — 30 
s\ \27r\s\J 



(AL/2vr) 2 - \s\ 2 



1*1 



Numerically we have for these three cases 



0.0618A-0.180L" 1 (cube) 
J 3 (0) = { 0.0573A - 0.314L- 1 (cilinder) , (31) 
0.0507A - 0.310L- 1 (sphere) 



4 In general, correspondence between £i£> and the minimal Lyapunov exponent is not so straightforward 
!]; in the present paper we ignore such complications. 
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Figure 3: Comparison of the theoretical scaling curve for d = 2 with numerical results by 
MacKinnon - Kramer [T8| Fig. 2, a] (a) and Schreiber- Ottomeier [191 Fig- 4] (b). 
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i. e. value of b\ is not universal but depends on the way of cut-off. The change of this 
constant allows to make the scaling curve more symmetric, or less symmetric; it was chosen 
from the best agreement!!, though its variation in the interval (—0.3, 0) does not affect the 
results significantly. As in the 2D case, the absolute scales for £ and £i£> are not fixed by 
the theory. 

Using the asymptotic behavior of H(z) 

f l/2z, z<l 
H(z)=l -A(z-z*), z^z* , (32) 
[ -c d z d - 2 , z > 1 

one has in variables y = £m/L and x = £/L 

2c d /x d ~ 2 , y > 1 

y*±B/x d ~ 2 , y^y* , (33) 

x , |/<Cl 

where 2* and y* = 1/z* are values of z and y in the critical point. The same relation, 
considered in variables y and 1/x, determines the L-dependence of the scaling parameter 
(Fig. 1) giving two universal curves for r > and r < to which all other curves are 
reduced by the the scale transformation: 

~ rL d ~ 2 , y > 1 

y = ^ = { y* + const ■ rL d ~ 2 , y y* . (34) 

In Fig. 4, the obtained scaling curves are compared with the early results by MacKinnon 
- Kramer [18] and the more precise results by Markos |IJ . In the former case the agreement 
is satisfactory, in the latter case there is discrepancy on the level of 2 - 3 standard deviations. 
However, one should have in mind how scaling curves are constructed: the L-dependences 
for different r are "measured" in the interval (L min , L max ), and then they are fitted to 
each other by a change of the scale (Fig. 5). The full scaling curve is never present in one 
experiment, and only separate fragments of it are measured. It is clear from Fig. 4,b that the 
change of the scale along the horizontal axis (reducing to a parallel shift in the logarithmic 
coordinates) allows to obtain satisfactory fits for the left, right or middle portion of the 
curve. It looks, there are no serious contradictions with the Vollhardt and Wolfle theory 
on the level of raw data. 



4. DISCUSSION OF THE SITUATION AT D = 3 

The interesting question arises: if the Vollhardt and Wolfle theory describes the raw 
data successfully, then why all numerical experiments give v > 1 for d — 3 ? 



' The fitting was made by hand (using several reference points) and probably is not optimal. 
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3D case 




Figure 4: Comparison of the theoretical scaling curves for d = 3 with numerical results 
by MacKinnon - Kramer [T8"l Fig. 2,b] (a) and Markos [TJ Fig. 53, right] (b). The values 
b\ = —0.240 and b\ = —0.0718 were used in the former and the latter case correspondingly. 
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Figure 5: Construction of scaling curves. 

The history of this question goes back to two papers [T7j and [IS] by MacKinnon and 
Kramer, based on the same array of the data. The first of them gives the result 

i/ = 1.2±0.3, (35) 

compatible with the value v = 1; the second paper confirms this result for a certain fitting 
procedure, but reports the "more precise" result 

v = 1.50 ±0.05, (36) 

corresponding to the most extremal of present-day values. The first result is based on the 
analysis of the scaling curve, whose compatibility with the Vollhardt and Woffle theory is 
clear from Fig. 4, a and is confirmed by the authors themselves. Further, they indicate that 
scaling is not satisfactory in the small vicinity of the critical point, and this vicinity was 
discarded in their treatment. In fact, such situation is natural, because the small vicinity 
of the transition is strongly affected by scaling corrections (see Eq. 23); the latter are small 
in magnitude but should be compared with the small value of r. However, the authors of 
[TTt IT8] estimated this situation as internally inconsistent and suggested another treatment 
procedure, which is specially based on the analysis of that small region where scaling is 
absent. Already at this stage it is possible to understand that the latter procedure is not 
satisfactory. 
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Indeed, using the systems of the restricted size L, one can work straightforwardly only 
in the regime £ < L, since in another case the correlation functions are strongly affected by 
finiteness of the system. The use of finite-size scaling allows "to jump above the head" and 
to advance in the region £^L; however, it is possible only if (a) scaling exists theoretically, 
and (b) it is observed empirically. If any of two conditions is violated, no such advancement 
is possible, and one is unable to obtain any experimental information on the large £ region; 
any manipulations in this region become irrelevant. This conclusion is valid in respect of 
the result (36), since absence of scaling is admitted by the authors. The same conclusion 
follows from the common sense: if value v — 1 is compatible with the scaling curve, then all 
the more it is compatible with the raw data (see the end of Sec. 3). However, this value is 
rejected by the result (36), and hence the latter should be qualified as essentially incorrect. 

The indicated tendencies was continued in other papers. The treatment based on the 
scaling curves gave rather conservative estimates^, not very different from (35). The results 
close to (36) were stabilized only when the control of scaling ceased to be imperative and 
the analysis of small vicinity of the critical point was generally accepted. 

The latter procedure is based on representation of (4) in the form 

i -f = F (^pj = F (tL^) *y* + AtL^ + . . . (37) 

i. e. the regular expansion in r is used, motivated by the absence of phase transitions in 
quasi- ID systems; then the derivative over r behaves as L x l v and gives the exponent v 
straightforwardly. Such treatment is correct if the scaling relation (4) is exact. However, 
it is not exact: linearization of (23) gives 



a 2 

T ~ c jr 

SID 



+ o(j) . (38) 



Differentating over r and excluding (£id)' t from the right-hand side in the iterative manner, 
one has 



/ 



(i^j = A L d ~ 2 + AiL d ~ 6 . (39) 

Producing subsequent iterations and taking into account further corrections to scaling, one 
have the following structure of the result 



L 



y* = t {AqL 1 ^ + A 1 L UJl + A 2 L^ + ...} + BJT^ + B 2 L^ + . . . , (40) 



which can be obtained from the general considerations based on the Wilson renormalization 
group [2]. 



6 As was discussed in [35], the estimate of v depends on the fragment of the scaling curve, which is used 
for fitting. 
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In three dimensions, the main scaling correction in (39) reduces to constant, and hence 

^- y * = Ar(L + L ) , (41) 

where the terms dissapearing at L — > oo are neglected. It is clear from Fig. 6, a, that 
numerical data by Markos [1] are excellently fitted by (41). The author himself interpreted 
them in accordance with (37) and also had the good fitting (Fig. 6,b). 

Such ambiguity of interpretation has a general character. If the combination A\L^ X + 
A 2 L^ 2 can be linearized in the log-log coordinates with the average slope (/?i + /3 2 )/2 and 
the accuracy e, then variation f5\ — > f3\ + 5, (3 2 — > fa — 8 preserves linearity on the same 
level of accuracy, till \8\ ^ — fa\/2. If several terms are retained in (40), the situation 
becomes not controllable at all: non-linear fitting with minimization of x 2 reveals the huge 
number of minima, and the deepest of them is not necessary correct; a vicinity of any 
minimum is acceptable if it is satisfied to the \ 2 criterion. Analysis of all such minima is 
impossible, and there is no honest procedure to deal with such situation^. In conclusion, 
the conventional treatment is heavily based on the assumption, that only the main term 
in (40) is essential; the problem of fitting becomes hopeless, if additional terms are not 
negligible. 

In the framework of the Vollhardt and Woffle theory, we have a completely consistent 
picture. The quantity Lq violates scaling and empirically has rather large value, Lo ~ 5 
(in lattice units). The good scaling is possible only for L ^> 5, and even for the largest 
systems (L = 20 -j- 30) deviations of scaling are described by the parameter Lq/L ~ 0.2; 
so discrepancies in Fig. 4,b should not be of any surprise. The theoretical value of L is of 
the order A -1 with the coefficient depending on the way of cut-off; it is essential that L is 
positive and limited from below by the atomic scale. 



5. SCALING FOR HIGHER DIMENSIONS 

5.1. Dimensions d > 4 

For d > 4, the sum ^{m) is divergent at the upper limit and the cut-off parameter A 
cannot be considered as infinite. For the accurate trasformation, we introduce the scale Ai 
such as 

m < Ai < A (42) 

and divide summation in ^{m) into two regions |g_i_| < Ai and \q±\ > Ai. In the first region 
we use that \q±\ <C A, so as 

J«( TO ) = -m 2 — ^ ]T 1 (43) 

2L q±jL0 \q±\ \Jm 2 + ql (|g±| + \jm 2 + ql) 

|?x|<Ai 



7 These questions were discussed [2] in relation with the paper [29]. Nevertheless, this paper is continued 
to be cited [T] as a prominent achivement. 
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Fitting by c(L+L ) 
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Figure 6: Numerical data by Markos for Z\ = 2L/£id in the small vicinity of the critical 
point [1, Fig. 53, left] and their fitting by dependences (L + L ) (a) and L - 63 (b). 
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and 

— m 



4V) 



-m 2 



— - — — + 0[m )\, m>L 
4(d ~ 4) . } (44) 

Ad-4 S J l^V 

^4i)- + 0(L 4 -")|, mS£ - 

i. e. the result is obtained analytically (in the main approximation) for the arbitrary relation 
between m and L^ 1 . Indeed, for m ^> L^ 1 the sum is estimated by the integral, which 
is converging at the lower limit already for m = 0; so a finiteness of m gives only small 
corrections. In the case m < L" 1 , the main effect from a finiteness of L is related with the 
absence of the term qj_ = 0, which can be estimated as restriction \q± \ <^L _1 in the integral 
approximation. 

In the region \q±\ > Ai we make use of condition \q±\ ^> m and produce expansions 
in m/\q±\; after separation the factor m 2 we can set m = in the sum and estimate it by 
transformation to the integral 

jy- a d— 4 

t(2)/ \ 2- fY d-l YV l 2\d-4 / AA \ 

In (m) = m — — cm A 44 

v ' 4(d-4) 

where c depends on the way of cut-off; dependence on Ai dissapears in the sum I^ + 1^ ■ 
The results for ii(m) and ^(0) are the same as in Sec. 3. The self-consistency equation 
takes the form 

rA d - 2 = J— JL _ cm 2 A d - 4 (45) 
Substituting r ~ £~ 2 and introducing variables 

^ = — UJ ' x = zUJ ' (46) 

we obtain the scaling relation in the analytical form 

4=»"7 (47) 

where all coefficients are made equil to unity by redefinition of the scales for £m and £. 
Relations (46 ,47) contain the atomic scale a, as was expected from non-renormalizability 
of theory (Sec. 1). 

According to (46,47), the role of the scaling parameter is played by the quantity y 
instead of £id/L; the Independence of y is analogous to Fig. 1, i. e. all curves corresponding 
to t > and r < can be reduced to two universal ones by the scale transformation. The 
transition point corresponds to y — 1, so as 

and the critical point cannot be fixed by the condition £id/L = const. 
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5.2. Four-dimensional case 



In the case d = 4 we have analogously 



h{m) 



-C4 m 2 In — + O (1) , m > L' 1 

m 

-c 4 m 2 ln(AL) + 0(1), m < 



(49) 



i.e. two results differ by ln(mL), which reduces to the double- logarithmic quantity in the 
actual region (see below). Neglecting such quantities, we can obtain scaling also for d — 4. 
The self-consistency equation has a form 



rA 2 



1 



2mL 3 



and after changing to £ and £id 



c 4 



±P ln(f/a) = ^ - ^ ln(£i D /a) 



21 A 

c 4 m in — 

m 



c 4 



(50) 



(51) 



The scaling relation (47) is obtained in variables 



IHL/a)} 



-1/3 



X = 



[HL/a)\ 1/e 



[ln(£/«)] 1/2 



(52) 



i.e. the scaling parameter y is logarithmically modified in comparison with £id/L an d 
should be considered as a function of "the modified length" = L [h^L/a)]^ 6 ; then 

all dependences become analogous to Fig. 1 and a change of the scale for /x(L) allows to 
reduce them to two universal curves for r > and r < 0. The critical point corresponds 
to y — const, so as 



ID 



T = 



(53) 



i.e. parameter £id/L grows logarithmically in the transition point. 

5.3. Modified scaling for d = 4 — e 

From the methodical point of view, it is interesting to derive the modified scaling for 
d = 4 — e; in this case, the sum ^ijn) converges formally at the upper limit, but this 
convergence is slow and a finiteness of A gives the essential effect. Analogously to (49) we 
have 

2 m - £ -A- £ 

, m> L 



him) 



-c 4 m 



,L e - A" 



-c 4 m 



m<L- 1 



(54) 
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and the scaling relation (47) is obtained in variables 



;id 



1 - (L/a) 



1/3 



x 



e^/a) [1 - (L/a)-] 1/6 



U/ a y - 1] 1/2 (L/ a y-^ 



(55) 



Once again we have the modified scaling parameter y and "the modified length" [i(L) 



-1/6 



in terms of which the dependences of Fig. 1 are recovered. In the 



critical point we have y = 1 and 



L 



1 - (L/a) 



1/3 



[\n(L/a)f 3 
(l/e) V3 , 



ln(L/a) < 1/e 
ln(L/a) > 1/e 



(56) 



i. e. parameter £ 1D / ' L grows logarithmically till the large length scale L ~ aexp{const/e}, 
and then saturates at the constant value. Such kind of scaling is useful as an alternative 
treatment for d = 3, in order to investigate the systematic errors related with the possible 
existence of the large length scale. In this case the parameter a in fact corresponds to Lq 
and can be essentially different from the lattice constant; it should be adjusted from the 
condition of the best quality of scaling. There is no need to be bound by Eq. 47, which 
is valid for small e; it is more reasonable to determine the relation y = F(x) empirically. 
As for expressions (55), their extrapolation to e = 1 does not present any problem, since 
for L, £ 3> a the modified scaling safely reduces to the usual one (see Eq. 4). In fact, it is 
identical to (4) if no large scale L$ is present. However, in the presence of the large length 
L such scaling is more adequate than (4). 



6. CONCLUSION 

The above analysis allows to conclude that the Vollhardt-Wolfle theory does not have 
essential contradictions with numerical results on the level of raw data. The different 
critical behavior reported usually in numerical papers originates from the fact, that some 
time ago the pure "experimental" approach to the problem was rejected and replaced by 
phenomenological analysis, which is practically hopeless in the corresponding region. In 
particular, dependence L + L with L > is interpreted as L x l v with v > 1. 

We have restricted our discussion by the widespread variant of finite-size scaling, based 
on application of auxiliary quasi- ID systems. Apart it, another algorithms are used, based 
on the level statistics [20J, the conductance distribution, the mean conductance, etc. pQ. 
The scaling curves calculated above are not universal and cannot be used for comparison 
with such results. The scaling for higher dimensions is also not universal: for example, 
another behavior in the critical point is expected for the Thouless parameter [40]. The 
only exclusion is the result (40), which remains unchanged in all cases. Indeed, this result 
can be obtained from the general arguments based on the Wilson renormalization group 
[2]; the exponents Ui, w 2 , yi, y 2 , ... are determined by scaling dimensions of irrelevant 
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Figure 7: Fitting by dependence (L + Lo) (dashed line) for numerical data, based on the level 
statistics: (a) Data by Zharekeshev and Kramer [27] • The points correspond to the average deriva- 
tives of the scaling parameter A (arbitrary units), determined from Fig. 4 of [27] in the interval 
16 < W < 17. A statistical error related with each point can be estimated very conservatively 
(see Table in [2]) due to the irregular character of curves given in the indicated figure; uncertainty 
allowed by the authors themselves corresponds to the gap between dependences L 0,80 and L 0,65 , 
determining the upper and lower bound of the result for the critical exponent, v = 1.40 ± 0.15. 
(b) Data obtained by Schreiber' group [33]; the points correspond to the derivative of the scaling 
parameter a (arbitrary units) determined by the slope of solid lines in the inset of Fig. 3 in [43J; 
their uncertainty is obtained by variation of the slope allowed by the size of experimental points. 
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parameters and hence are universal. Correspondingly, the result (41) is unchanged, which 
explain the origin of the effective values v > 1 (Fig. 7). 

Fig.7,a can be considered as an etalon illustration, corresponding to most of numerical 
papers. Indeed, there is an overall consensus that data for L < 5 fall out of the scaling 
picture and should be discarded; large systems with L > 30 are practically never used; the 
error corridor between dependences L 0,80 and L - 65 corresponds to the typical accuracy of 
numerical papers. Fig.7,b illustrates one of the rare papers treating the systems of the 
record size [43] . Finally, Fig. 8 shows the rare example of high-precision data |44j . 

Our final remark is as follows. Even if subsequent investigations reveal that the Voll- 
hardt and Wolfle theory is not exact, nevertheless no confidence can be given to the present- 
day estimates of the exponent v [L7] [3T|- Fig. 6 clearly demonstrates that values i/wl.6 
and v = 1 are equally compatible with the raw data, and hence the treatment procedure 
is extremely ambiguous. 
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